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Abstract In this article, we are interested in the large time behavior of solutions of 
the Dirichlet problem for subquadratic viscous Hamilton-Jacobi Equations. In the su- 
perquadratic case, the third author has proved that these solutions can have only two 
different behaviors: either the solution of the evolution equation converges to the solution 
of the associated stationary generalized Dirichlet problem (provided that it exists) or it 
behaves like — ct + ip(x) where c > is a constant, often called the "ergodic constant" and <p 
is a solution of the so-called "ergodic problem" . In the present subquadratic case, we show 
that the situation is slightly more complicated: if the gradient-growth in the equation is like 
\Du\ m with m > 3/2, then analogous results hold as in the superquadratic case, at least if 
c > 0. But, on the contrary, if m < 3/2 or c = 0, then another different behavior appears 
since u(x, t) + ct can be unbounded from below where u is the solution of the subquadratic 
viscous Hamilton-Jacobi Equations. 

Key-words : Viscous Hamilton-Jacobi Equations, large time behavior, subquadratic case, 
Dirichlet problem, ergodic problem, viscosity solutions. 

AMS subject classifications : 35K55, 35B40, 49L25 



*Laboratoire de Mathematiques et Physique Theorique (UMR 6083). Federation Denis Poisson 
(FR 2964) Universite de Tours. Faculte des Sciences et Techniques, Pare de Grandmont, 37200 
Tours, France. E-mail address: barles@lmpt.univ-tours.fr. 

^Dipartimento di Matematica, Universita di Roma Tor Vergata, Via della Ricerca Scientifica 1, 
00133 Roma, Italia. E-mail address: porretta@mat.uniroma2.it. 

■^Departement de Mathematiques, Faculte des Sciences, Universite de Yaounde I, BP: 812, 
Yaounde, Cameroun. E-mail address: tabet@lmpt.univ-tours.fr. 



1 



Contents 



1 Introduction 

2 Preliminary Results on the Stationary Ergodic and Dirichlet Prob- 
lems 

3 Asymptotic Behavior for the Parabolic Problem |ll 

3.1 Convergence to the solution of the Stationary Dirichlet Problem ... Ill 

u(Xj t) I — 

3.2 Convergence of — 12 

t 

3.3 Convergence to the Stationary Ergodic Problem when § < m < 2 and 



c> 14 



4 The Non-Convergence Cases: c = or c > and 1 < m < 3/2 [20 



5 Appendix [25 



1 Introduction 

We are interested in this work in the behavior, when t — > +oo, of the solution of the 
following initial-boundary value problem 

u t -Au+\Du\ m = f(x) inftx(0,+oo) (1.1) 
u(x, 0) = uo(x) on Q (1.2) 
u(x,t) = g(x) on dQ x (0, +oo) (1.3) 

where Q is a C 2 - bounded and connected subset of ~R N , 1 < m < 2 and f,u ,g are 
real-valued continuous functions defined respectively on Q and <9fi. The boundary 
and initial data satisfy the following compatibility condition 

uq(x) = g(x) for all x G dVl. (1.4) 

It is standard to show that this problem has a unique solution u : Vt x [0, +oo) — * R 
and, as long as one does not need regularity properties, proofs of this fact are easy 
by using viscosity solutions theory (see Barles and Da Lio [1] and references therein). 

In the superquadratic case (m > 2) the study of the asymptotic behavior has been 
done by the third author in [16] where it is shown that the solution u can have only 
two different behaviors whether the equation 

- Av + \Dv\ m = f{x) in ft, (1.5) 
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has bounded subsolutions or not. If (II. 5p has a bounded subsolution, then there exists 
a solution of the stationary Dirichlet problem, i.e. of (11.51) together with the 
generalized Dirichlet boundary condition 

v(x) = g(x) on dQ , (1.6) 

and u(x,t) — > Uoo(x) uniformly on Q. This is the most expected behavior. On the 
other hand, it can happen that (II. 5p has no bounded subsolution and, in this case, 
one has to introduce the ergodic problem with state constraint boundary conditions, 
namely (the reader has to keep in mind that we are here in the super quadratic case) 

f -A<p+\D<p\ m = f(x) + c inn, 

\-A(p+\D(p\ m > f(x) + c onffi. 1 ' 

We recall that, in this type of problems, both the solution (p and the constant c 
(the ergodic constant) are unknown. The existence and uniqueness of solutions (c, <p) 
for ( II. 7j) is studied in Lasry and Lions [H] (see also [H] for a viscosity solutions 
approach): the constant c is indeed unique while the solution p is continuous up to 
the boundary and unique up to an additive constant. 

Concluding the reference when m > 2, it is proved in [16] that, if 1I1.5P has no 
bounded subsolution, then c > and the function u(x, t) + ct converges uniformly on 
O to a solution p of 11 1.7ft when t — > +oo. In fact, in the superquadratic case, even if 
the boundary condition reads u(x,t) + ct = g(x) + ct on dfl with g(x) + ct — > +oo 
as t —>■ +oo, there is a loss of boundary condition (cf. [1]) and u(x,t) + ct remains 
bounded on Q. 

The first key difference in the subquadratic case is that there is no loss of boundary 
conditions and u(x,t) + ct is actually equal to g(x) + ct on the boundary. Formally 
this forces the limit of u(x, t) + ct to tend to +oo on the boundary and, actually, the 
analogue of 1I1.7P is 

f -A<p+\D<p\ m = f(x) + c inn, 

\ p(x) — ► +oo when x — > dn . 

This problem was also studied in [TJ] where it is proved that there exists a unique 
constant c such that 1I1.8P has a solution and, as in the superquadratic case, the 
solution is unique up to an additive constant. 

Coming back to the asymptotic behavior of u, it can be thought, at first glance, 
that it is essentially the same as in the superquadratic case, with problem 111.71) being 
replaced by problem fll.8p . Surprisingly this is not true for any 1 < m < 2 nor for 
any c. Part of the explanations concerning c is that the equation 

-Aip+ \Dip\ m = f{x) + c inn, 
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has subsolutions which are bounded on Q for any c > c in the superquadratic case 
but only for c > c in the subquadratic one. 

In this article we prove the following: first, if the stationary Dirichlet problem 
fll.5p - fll.6l) has a solution, then u converges uniformly in Q to this solution. Oth- 
erwise we show that necessarily c > and u can have different types of behavior. 
Observe that a curious feature in the subquadratic case (which does not occur in 
the superquadratic case) is that c can be equal to even if the stationary Dirichlet 
problem has no solution; this can be seen as a corollary of the above remark. The 
different behaviors of u can be described as follows. 

(i) If 3/2 < m < 2 and c > then u(x, t) + ct converges locally uniformly in Q to 
ip where <p is a solution of (11.81) . 

u(x t) 

(ii) Ifl<m<3/2 and c > then — j > — c locally uniformly in f2 but it can 

happen that u(x, t) + ct —* —00 in Q. 

u( x t) 

(iii) If c = 0, then — - — — > locally uniformly in Q but it can happen that 
u(x, t) — > —00 in Q. 

The behaviors (ii) and (iii) are striking differences with the superquadratic case. 
This is related to the blow-up rate of u(x,t) which we estimate in Theorem 13.21 and 
which is influenced by the behavior of p (solution of fll.Sp ) near the boundary dfl 
(changing according to the values of m) and eventually by the case that c = 0. We 
prove the optimality of such estimates in star-shaped domains in cases when / + c < 
on Q: in such situations, we show that it can actually happen that u(x, t)+ct — > —00. 

This article is organized as follows. In Section [2j we recall the main results about 
the problem (jl.8P ; these results are mainly borrowed from Lasry and Lions [14J and 
Porretta and Veron [15] where the precise behavior of <p near dfl is described. Sec- 
tion EH] is devoted to the easy case when (11.5p - (11.6H has a solution; in such situation 
u(-,t) is bounded and converges to the (unique) solution of the stationary problem. 
We recall that the existence of a stationary solution corresponds, for the ergodic con- 
stant c of (11.81) . to the case c < (see Proposition 12. 2p . In Section 13721 we consider 
the case c > and we state the estimates on the blow-up rate of u proving, as a 
consequence, that "^'^ — > — c locally uniformly in f2, in the whole range 1 < m < 2. 
In the Section I3T31 we prove the result (i) above. Finally, in Section 0] we show that a 
similar result cannot hold in general in the range m > | or if c = proving that, in 
these cases, we have u + ct — ► — 00 at least in some circumstances. We leave to the 
Appendix the construction of sub and super solutions, for the problem fll.ip - fll.2p - 
fll.3p . which we use to obtain the crucial estimates on the blow-up rate of u when 
c > 0. 
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Throughout this paper, we assume that Q is a domain of class C 2 . We denote by 
d the signed-distance to dQ, which is positive in Q, i.e. d(x) := inf{|x — y\ : y £ dQ} 
if x E Q and negative in the complementary of Q in M. N . As a consequence of the 
regularity of dfl, d is a C 2 - function in a neighborhood W of <9f2. We also denote by 
1/ the C 1 - function defined by u(x) = —Dd(x) in W; if x e 9f2, then v(x) is just the 
unit outward normal vector to dVt at x. 

2 Preliminary Results on the Stationary Ergodic 
and Dirichlet Problems 

In this section, we deal with properties of the pair (c, ip) solution of (11.81) . We start by 
recalling the main results concerning (c, ip) which were proved in [TJ] and \J5\. Next, 
we provide some useful properties of c giving a relationship between the stationary 
ergodic problem (11.81) and the stationary Dirichlet problem fll.5p - fll.6l) . We end the 
section with an example to illustrate the possibly unsolvability of (|1.5|) - ()1.6p . 
We start with 

Theorem 2.1 (On the Stationary Ergodic Problem) 

Assume that 1 < m < 2 and that f G C(Q). There exists a unique constant cGl 
such that the problem U.8\) has a solution (p e Wf^iVl) for every p > 1. This solution 
is unique up to an additive constant and satisfies the following properties 

(i) If 1 < m < 2, then 

2 — m f rn — 1 ) m ~ 1 

<p(x) = C*d-—i(x)(l + o(l)) as d(x)^0 , with C* = '—— . (2.1) 

(2 — m) 

(ii) Ifm — 2, then 

<p(x) = \]ogd(x)\{l + o(l)) asd(x)-+Q. (2.2) 
(Hi) There exists a constant A > depending only on m, Q and f such that 

\D<p(x)\ < Ad~^(x) in Q. (2.3) 

(iv) We have 

lim d~^ = i(x)D(p(x) = c m v(x) (2.4) 

where c m = (m — l) - ™- 1 . 
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As far as the proof of Theorem 12.11 is concerned, we refer the reader to the proofs of 
[TH Theorems 1.1, IV. 1 and VI. 1] for (i), (ii), (iii) and to [ISl Theorem 1.1-B] for (iv). 

We continue with further estimates on the solutions of (jl.8p . which can be deduced 
by the previous ones. 

Lemma 2.1 Assume 1 < m < 2, f G W l '°°(Q) and let if be a solution of U.8\) given 
by Theorem^ Then if G C 2 f c (Vt) for all (3 G (0, 1) and 

(i) There exists a constant K > such that, for every x G Q, we have 

\D 2 <p(x)\ < K d- a - 2 (x) with a = . (2.5) 

m — 1 

(ii) There exists K\ > and a > such that, for every x G Q with < d(x) < a 

\D\{x)\<K x \D^{x)\ m (2.6) 

and 

\<p(x)\ < K x \Dp(x)\ 2 - m ifl<m<2, 

(2.7) 

|^(o;)| < \ log(\D<p(x)\)\ + K 1 ifm = 2. 

Proof of Lemma 12.11 We first remark that the C^'f-regularity of the solutions 
of (11.81) comes from the additional assumption on / (/ G iy 1,00 (f2)) and a standard 
bootstrap argument : indeed, by Theorem 12 .11 ip is in Wf^{Vt) for every p > 1, hence 
in C]£{Vt) for any (3 G (0, 1), and therefore a standard regularity result implies that 
if G Cff c {Q) for all p G (0, 1) since \Dp\ m and / are in C°f (ft). 

In order to prove (i), let Xq G fi, and r = We introduce the function 

ip(x) = r a (p(xo + rx) and compute 

D^(x) = r a+1 Dip(xQ + rx) and D 2 ip(x) = r a+2 D 2 ip(x + rx). 

Since m(a + 1) = a + 2, it is easy to see that if) is a viscosity solution of 

-Aip + \Di)\ m = r a+2 (f(x + rx) + c) for any x G 5(0, 1) . 

Since r a+2 f and \Df\r a+2 are in L°°(f2), we can use the interior estimates available 
in [TU Theorem A.l] and obtain 

IW IL- ( b( 0) i)) < ^ 
and then, by elliptic regularity (see e.g. [13]), 

11^11^(5(0,1)) < ^' 
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which yields 

\D 2 ip(x + rx)\ < K'r- a - 2 for all x E 5(0, 1/2). 

By taking x = and remembering the definition of r, we obtain (12.51) with K := 
2 a+2 K'. The estimate (12.61) is a consequence of (12.41) which implies 

\Dip(x)\ > Kd~^(x) (2.8) 
for some constant K > and for x in a suitable neighborhood of dfl. Since — ^ = sdl 

° m— 1 m ' 

combining (12.81) with (12. 5p we obtain 

|DV(x)| < K d~ a - 2 (x) = K d-^(x) < K (^\Dv(x) 

hence i \2.6h for any K\ > K /K m . As far as (12. 7h is concerned, we use (12.1 ft and (12. 8ft 
if m < 2 and (12.21) and (12.81) if m = 2 obtaining that (12.71) holds true for some K\ > 0. 

□ 

We continue by showing monotonicity and stability properties of the ergodic con- 
stant c with respect to the domain. 

Proposition 2.1 

(i) Let Q' be an open bounded subset ofH N such that Q C Q'. Let cq and cqi be the 
ergodic constants associated to U.8\) in Q and Q' respectively. Then we have 

cn < en*. (2.9) 

(ii) Let Q' be an open bounded subset ofR N such that O CC O'. The respective ergodic 
constants cn and Cqi of Q and Q' satisfy 

c n < c Q/ . (2.10) 

(Hi) Moreover, the ergodic constant c depends continuously on Q. Otherwise said, for 
< i] < 1, if Cjj is the ergodic constant in JZHP set inVt + £> (0,^)0 then 

c v i c n as r) i (2.11) 

where B(0,r]) is a ball of radius r\. 



x ft + 5(0, t?) := {x G R N : d(x) > -n}. 
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Proof of Proposition |2~TTI (i) Let (cq, (pn) and (cqi, </pq') be the pair of solutions of 
the ergodic problem (11.81) in fl and fl' respectively. From what we obtained above, 
the constants cn and cw are unique whereas the functions (pn and (pn> are unique up 
to a constant. We study max^-(yU<^/ — (p&) for some /i G (0, 1) close to 1. Now, we 
need to show that the maximum is achieved inside fl. 

Observe that, since fl C fl', it could happen that dfl fl dfl' ^ 0, meaning that the 
two domains fl and fl' touch at some points. Due to the behavior of (p^ and ifn near 
dfl' and dfl respectively, we first deal with points on dfl which do not belong to dfl' 
and next we treat common points of dfl and dfl'. 

We pick any xq G dfl such that xq ^ dfl'; in this case, since <pw is bounded in 
H \ (dfl n an'), it follows that 

lim (/iy?fy - <^n)(x) = //</?fi'( x o) - lim </?n(z) = -oo. 

On the other hand, if xq G dfl fl <9fi', using the asymptotic behavior (12.11) and that, 
if daw denotes the distance to dfl', ddw(x) > d(x) we have, when 1 < m < 2 

2 — - ui 

(H<fn> - <pn)(x) < C*(/i - l)d~^(x)(l + o(l)) asx^xo, 

hence (fnpn' — fn)^) — > — oo in this case too. When m = 2, the same conclusion 
holds by the use of (12. 2p . Therefore, in any case it follows that (fupn> — ipn) has a 
maximum point x M G fl. Going back to the equations solved by fnpw and (p^, we 
obtain: 

- A(/^ n 0(x M ) + |£(my)Mr = -^A^(x M ) +/i m |D^(^)| m 

•^m / 1 / 

= fifM+fica, (2.12) 

and 

- A^(x M ) + |ZV n (z M )| m = f(xj + en- (2.13) 

By subtracting (12. 13[) from (12.121) and using the following properties: D(fi(pQi)(x^) = 
Dip n (x IJ ) and A(fap n > - (pn)(x^) < 0, one gets: c Q < fic Q > + (/i - l)f(x^) and §ZM) 
follows by sending jtt to 1. 

(ii) Let fl' be such that fl CC fl'; from (12.91) it follows that < cq'. Moreover, we 
would like to show that for this case, we have: cq < cqi; to do so, we assume the 
contrary by setting cn = cjy. Using the boundedness of (pn> on fl and behavior of <pn 
near dfl, it is easy to see that (fa — (pn — > — oo on dfl, meaning that (pa —<pn achieves 
its global maximum on fl at some x inside fl. On the other hand, the convexity of 
p i — ► \p\ m yields 

\D?n>\ m > \D<p n \ m + m\Dp n \ m ~ 2 Dp n ■ (D( m , - <p n )). 
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With this argument and using the local bounds on |£Vn| in (12.31) . we have that 
tpu> - tpn solves 

- A(cpw - (pa) — C(x)\D(cpQ f — (po)\ < in n (2.14) 

for some C(x) > which is bounded in any compact subset. Applying the Strong 
Maximum Principle (see Lemma 2.1]) we find that ipni — (fQ is constant in Q 
which clearly leads to a contradiction since pn> — (fn blows up on the boundary dfl, 
whence (12.101) holds. 

(iii) Now, we turn to the proof of the continuous dependence of c and p in Q, namely 
(12.111) . Let < t] < i]' < 1, since + 5(0, rj) is a strict subset of Vt + 5(0, 77'), we find 
from (12.1 Op that c r] < c v >, meaning that the sequence (c r) ) <r ) <i decreases as r\ goes to 
0. On the other hand, knowing that Q C Q + -8(0, rj) for any < 77 < 1, again by 
(12. 9p . we obtain cn < c v . Therefore, (c v )o <rj< i is convergent in R and { cq = cq as 
1] — ► + by using the arguments on stability of (11.81) and the uniqueness of the ergodic 
constant. □ 

Hereafter, to stress on the dependence of the initial boundary value problem ( 11.11) - 
( jl.2p - lll.3p on /, m and g, we denote it by E(Q, f, g, « ) where "5" stands for Evolu- 
tion. Likewise, we denote by S(Q, f, g), the boundary value problem fll.5p ~ ljl.6p where 
"5"' stands for Stationary. 

Now, we link the problems (11 ,8p and S(Q, f, g) by pointing out that the existence 
of a solution for S(Q, f,g) depends on the ergodic constant c in (11. 8p . 

Proposition 2.2 Let c be the ergodic constant associated to U.8\) and let us denote 
by S the set of all A G M such that there exists a viscosity subsolution G C (Q) of 

- A0+ \D(j)\ m < f + X inQ. (2.15) 

Then 

c = inf{A:AG5} (2.16) 

and the infimum in A2.16}) is not attained. Moreover, a necessary and sufficient 
condition for S(Q, f, g) to have a viscosity solution is that c < 0. 

Remark 2.1 It is worth mentioning that by a viscosity solution of S(Q,f,g), we 
mean a C(Q)- function satisfying U.5\) in the viscosity sense and U.6\) pointwisely. 
Indeed, when 1 < m < 2, there is no loss of boundary conditions as specified in |^} 
Propositions 3.1 and 3.2], and S(Q, f,g) is a classical Dirichlet problem. 

Proof of Proposition 12.21 In order to prove (12.161) . we first remark that 5^0; 
indeed, it is easy to see that ||/||oo G S because <p = is a subsolution of 

-A^+|£vr</+imioo m a 
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On the other hand, we pick any A G S and denote by (p\ G C(TT), the function 
satisfying — A0 A + \D(j)\\ m < f + A in the viscosity sense. We argue as in the proof 
of Proposition 12.11 by studying max.^(<f)\ — tp), in order to reach the conclusion that 
c < A. Therefore, the right hand side of f 1 2 . X 6 [) . which we denote by A*, is well-defined 
and c < A*. 

Now, we assume that c < A*, by keeping the notations of Proposition l2.H -(iii). 
there exists some < r] <C 1 such that, for any < r\ < r] , we have G]c, 
But, if ip v solves 

-Aip v + \Dtp v \ m = f + c v mQ v :=n + B(0,T]) 

then <p v is in Wf£ for any p > 1 and therefore (p v G C(Q); thus (p v is a viscos- 
ity subsolution of the equation associated to in Q and c v G S. This clearly is a 
contradiction since c v < A* and we conclude that c = A* and then (12.161) holds. 

Moreover, the infimum in (12.161) is not attained: indeed, if (p c G C(f2) is a subso- 
lution of the c-equation (12.151) and if <p is a solution of the ergodic problem (11.81) . then 
the maxjy ((p c — (p) is achieved at some point of Q since ip blows up on dQ and applying 
the Strong Maximum Principle exactly as in the proof of Proposition 12. II (ii) . we find 
that tp c — <p is constant in Q, a contradiction since (p c G C(f2) and (p(x) — > +oo as 
x — > dVl. 

Now we turn to the proof of the second part of the result. If the generalized 
Dirichlet problem S(Q, f,g) has a bounded viscosity solution, then G S, and then 
the first part implies that necessarily c < 0. 

Conversely, we assume that c < and prove that S(Q, f, g) has a unique bounded 
viscosity solution. To do so, we are going to apply the Perron's method (cf. [IT] . [8] 
and [9]) and in order to do it, we have to build sub and supersolution for S(Q, f,g). 

For the subsolution, since c < 0, we find from (12.1 ip (see Proposition 12. ip . the 
existence of < r] Q < 1 such that c v < for all r\ < r] where c v is such that there exists 
a function ip v which is a viscosity solution of the ergodic problem (11.81) in Q + B(0, rj) 
for all rj < rj . 

Since g is bounded on dfl whereas (p v is bounded in Q, there exists K > such 
that ip v — K < g on dVt. Therefore we have, at the same time, a subsolution of 
S(Q, f, g) required in the Perron's method, but also a strict subsolution of S(Q, f, g) 
which implies a comparison result for S(Q, f,g) (See [T6l Theorem 2.3]). 

For the supersolution, it is easy to see that l(x) = \x — x | 2 + HgHoo + 1 is a superso- 
lution of S{Q, f, g) for some x G R N such that B(x , (||/|| 00 + 2A^) 1 / m ) nH = 0. Using 
the (strict) subsolution <p v — K built above, the existence of the solution therefore 
follows by combining the comparison result and the classical Perron's method. □ 

We end this section by giving an example showing that the stationary Dirichlet 
Problem S(fl, f, g) does not always have a solution. 
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Example 2.1 We fix some R > and consider the one- dimensional equation 



- v " + \r]'\ m = -C m in (-R,R) (2.17) 

with C > 0. If r\ solves (I2.17p . then after some easy change of variable and computa- 
tions, we find that 

_ i'w , 
1 [ c ds 

x. 



C m ~ l hm. \s\ m + 1 

c 11 



It follows that 



n'lQ) \s\ rn + 1 - /_„ S m + 1 



Therefore, letting x — > -R, we obtain : 

C" 1 " 1 < — (2.18) 

Since R > is fixed and m > 1, by choosing C > large enough, we find that the 
inequality (12.181) cannot hold and we conclude that the ordinary differential equation 
(12.171) is not solvable for large C. Therefore, (11.5p - (jl.6j) is not always solvable as 
specified above when we considered / := — C m < 0. □ 



3 Asymptotic Behavior for the Parabolic Problem 

This section is devoted to the description of the asymptotic behavior of the solution 
u of the initial boundary-value problem E(Q, f, g,uo). 

3.1 Convergence to the solution of the Stationary Dirichlet 
Problem 

We start with the case that u converges to the solution of the stationary Dirichlet 
problem. The main result of this subsection is the following. 

Theorem 3.1 Let 1 < m < 2, / e C(0), u G C(Q) and g e C((90). Assume that 
S(Q, f,g) has a unique viscosity solution which we denote by u^. Let u be the unique 
continuous viscosity solution of E(Q, f, g,uo). Then, as t —>■ +oo, 

u(x, t) — > u OQ (x) uniformly for all x G fl. (3.1) 
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Proof of Theorem 13.11 We first notice that u is uniformly bounded on Q x [0, +00). 
Indeed, since solves S(Q,f,g), it is straighforward that solves E(Q, f, g, w^) 
and by [HI Corollary 2.1], one gets 

|| u(x^ t) Uqq I) oo || Uq Uqq II 00 • 

Next, from the uniform boundedness of u obtained above, we use the half-relaxed 
limits method to say that the functions 

u(x) = limsup u(y,t) and u(x) = liminf u(y,t) 



y-,x 

t i- + 00 



t — V + OC 



are respectively subsolution and supersolution of 01.5l) -( fL6l) . By definition of the 
half limits, we have u < u on Q but given that the ergodic constant c is strictly 
negative (see Proposition 12. 2p . there exists a strict subsolution for S(Q,f,g) and we 
can therefore apply [T6J, Theorem 2.3] to obtain u < u on Q. It is worth noticing that 
since no loss of boundary condition could happen (see [H Propositions 3.1 and 3.2]), 
the Dirichlet condition (11.61) is understood in the classical sense: u < g on dVt and 
u > g on dfl. Thus, we obtain = u = u on fl, meaning that (13.11) holds. 

□ 

3.2 Convergence of ■ — 

& t 

In this subsection, we assume that c > 0. Indeed, due to Proposition 12.21 the case 
c < is already described by Theorem 13.11 The main result of this section is the 
following convergence of "^'^ to — c which always holds true, even if at different rates. 

Theorem 3.2 Let 1 < m < 2, f e W 1 ' 00 ^), u G C(Tt) and g e C(dQ). Let u be 
the unique continuous viscosity solution of E(Q, f, g,uo). Then we have 

—j ' ~^ ~ c locally uniformly onVt as t —>■ +00 . 

In particular, for any compact set K C £1 there exists a constant Mk such that, as 
t — > +00, 



(i) if c> then 



' W^f 1 + c\\ c{ K) < ^ when\<m<2 
^ + c\\ c{ K)<M K ^ whenm = \ (3.2) 



1 ^ iA + c\\c(k)<^t when Km <| 



t , , — 
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(ii) if c = then 



u(x,t) I 



c{K) < Mk when m = 2 



I r*(ac,i) ii ^ M K 

\—t-\\C(K) S 



when 1 < m < 2. 



(3.3) 



Theorem 13.21 is a direct consequence of the following estimate. For technical 
reasons, it will be convenient to consider the unique solution <p Q of (\1.8h satisfying 
min^(p = 0. In this way, as said above, any solution of (11.81) is described as 
if = ifQ + k for some constant k. 

Theorem 3.3 Let Q be a domain of class C 2 . Let f 6 W 1,OQ (Q). Then we have, as 
t — > +oo 

(i) If c> there exists a constant M > stzcn that 

{u + ct> 7(t)y>o(z - /i(t)n(x)) — M i/ | < m < 2, 

u + ct > 7(*)^ (z - ii(t)n(x)) - M logt i/m = |, (3.4) 
m + ct > 7(t)v9 (x - /i(t)n(x)) - Mt^=^ if 1 < m < |, 

where n(x) is a vector field such that n(x) ■ Dd(x) < and 7(t), //(£) are positive 
continuous functions such that j(t) is increasing and ^(t) | 1 ; //(£) zs decreasing and 
Hit) I as £ — ► oo. 

(ii) If c = £nere exists a constant M > snca £/ia£ 

w > 7(*)Vo(aJ - ~ Mlog£ ifm = 2, 

u > j(t)cp (x - fx(t)n(x)) - Mt 2 ~ m ifl<m<2, °"' 

where j(t) j 1 ; nit) { as £ — > +oo. 

Knowing that <po is nonnegative and solves (11.81) . it is easy to observe that (po + 
||wo||i°°(n) — ct is always a super-solution of (ll.ll) - (ll.2p - (ll.3l) . Therefore, by means of 
Strong Comparison Principle for E(Q, f, g,uo) (see [H Theorem 3.1] or [161 Theorem 
2.1]), we have for all 1 < m < 2, 

u(x, t) + ct < (po(x) + \\u \\L°°(n) in Q x [0, +oo). (3.6) 

Theorem 13.31 provides some estimates on u(x,t) + ct when c > which can be used 
to locally bound u + ct from below in order to complement (13. 6p . 

We refer the reader to the Appendix for the proof of Theorem 13.31 and turn to the 
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Proof of Theorem 13.21 It follows from (13.61) that, in any compact set K, we have 
« + c < Mjl f or 

some constant Mx- On the other hand the estimate from below 
varies according to the values of m and whether c = or c > 0. Using the estimates 
(13. 4p and (13.51) in Theorem 13.31 and taking into account that ipo is locally bounded, 
we immediately deduce (I3.2p and (13.31) and in particular that ^2-^1 — > — c locally 
uniformly as t — > +oo. 

□ 

In the case that c > and 1 < m < |, orifc = 0, the rates of convergence given 
above cannot in general be improved, as we will see later (Theorem 14.11) . 

3.3 Convergence to the Stationary Ergodic Problem when 

| < m < 2 and c > 

The goal of this section is to describe the asymptotic behavior of the solution u of 
the generalized initial boundary-value problem (ll.ip - (ll.2l) - (ll.3p in connection with 
the stationary ergodic problem (II. 8p . Our main result is the following: 

Theorem 3.4 (Convergence result) 

Let f G W 1 ' 00 ^), u G C(O) and g G C(dQ). Let c G R and y? e W^(n) be the 
unique viscosity solution of the ergodic problem U.8\) such that min^ (f = 0. Let u 
be the unique continuous viscosity solution of E(Q, f, g,uo). Assume that c > and 
| < m < 2, then we have 

u(x, t) + ct — > ipo(x) + C locally uniformly in fl, as t — >• +oo (3.7) 

for some constant C depending on Q, f, c, Uq, and g. 

Recalling that all solutions of problem (I1.8P only differ by addition of a constant, 
one can rephrase (13.71) saying that u(x,t) + ct converges to a solution of the ergodic 
problem (11.81) . We choose to represent all solutions as <po + C in order to emphasize 
that there is precisely one constant C, depending on the data, which determines the 
asymptotic limit of u + ct — tp^. 

We also notice that by combining (13.61) and (13. 4p . it follows that 
^y(t)ipo(x — n(t)n(x)) — M < u(x, t) + ct < (po(x) + \\uo\\L°°(n) in fl X [0, +oo). (3.8) 
If lim (u(x, t) + ct) exists, then by sending t — > +oo in (13. 8p . we obtain 

ipo(x) — M < lim (u(x,t) + ct) < (p (x) + ||«o|U°°(n) 
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From the blow-up behavior of ipo near the boundary, it obviously follows that 

lim (u(x, t) + ct) — > +00 as x — > <9fi (3.9) 

which gives the expected behavior near the boundary since there is no loss of boundary 
condition for all 1 < m < 2. 

The rest of this section is devoted to the proof of Theorem I3.4L 

Hereafter, for all (x,t) G Q x [0, +00), we set 

v(x,t) := u(x,t) + ct and iw(x, £) = v(x, t) — <fo( x )- 

The function v(-, ■ + t) solves E(Q, f + c, g + ct, v(-, t)) whereas ipo is a supersolution 
of E(fl, f + c,g + ct,ipo) since <^o solves (11.81) . From the comparison principle, one 
gets, for all x G Q and s > t > 0, 

max w(x, s) < max w(x, t). 

It follows that the function t t— > m(t) := max^^ w(x, t) is non-increasing. Moreover, 
as a by-product of (13.8p . m(t) is bounded. Therefore m(£) J, m as £ — > +00. 

Now, since we want to deal with bounded functions and clearly w(x, t) — > —00 
when x — > <9f2, we choose any constant X > |m| and set 

z(z, t) := sup[u>(x, t), — if]. 

We notice that, since m(t) is bounded from below by m, we still have 

m(t) := maxw(x,t) = max[sup(w(x, t), — K)] = maxz(x,t) for all t > 0, 

and from (13. 8ft . z is uniformly bounded on fix [0, +00). Moreover, since w(x, t) — > —00 
as x — > dQ for all t > 0, we have z = — on <9f2 x [0, +00) and we also remark that 
z is a viscosity solution of 

( 4> t (x,t)-A(f)(x,t)+n(x)-D(f)(x,t) <0 infix(0,+oo) , . 

\ 0(x,t) + fsT =0 on <9f2 x (0, +00) 1 j 

where 

= mlD^o^r" 2 ^^)- (3.11) 

Indeed, by the convexity of p 1— > |p| m , w is a subsolution of the above equation and 
we recall that the maximum of two subsolutions is a subsolution. 
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In order to have an equation with continuous coefficients, we introduce the oper- 
ator 

C<p(x, t) := d(x)(j> t (x, t) - d(x)A<j>(x, t) + d(x)H(x) ■ D<p(x, t) , (3.12) 

just obtained from the previous equation by multiplying by d(x). From (12.41) . it is 
easy to note that 

\D<p {x)\ m - 2 D(p (x) = -{m-l)- 1 ^^ + o(-^) as d(x) -> 0. (3.13) 



d(x) \d(x 

It is therefore obvious to see that 7i has a singularity on dfl whereas d(x)Ti,(x) can 
be extended as a continuous function on Q. Using this new operator, we have 

Cz < in x (0, +oo) , 

which replaces the subsolution property in (13.101) . 

To complete the proof of (13.71) . we first give the following local Holder continuity 
of the unique solution v of E(Q, f + c, g + ct, u ), with respect to its x and t variables. 



Proposition 3.1 Let 1 < m < 2 and v be the unique continuous viscosity solution 
of E(Q, f + c, g + ct,u ). Then 

(i) For all r] > and for all v G (0,1), we have v(-,t) G C^(Q) for all t > rj. 
Moreover, if Kg := {y G Q : d(y) > 5} for any 5 > 0, then, for all t > r), the 
C°' u -norm ofv(-,t) on K$ depends only on v, 5, r], \\f\\oo an d the L°°-norm of 
v on K 8 /2 x [|, +oo). 

(ii) For any x G Q, v(x, •) G C^ ' c 2 (0, +oo) for any r] > 0. Moreover the C°'%-norm 
ofv(x,-) in [rj, +oo) depends only on z/, r], \\f\\oo o,nd the L°°-norm of v on 
Kd(x)/2 x [77/2, +00). 

We postpone the proof of Proposition 13 . 14 - (i) to the Appendix and refer the reader to 
[31 Lemma 9.1] for the proof of Proposition 13. II (ii) . 

In order to prove (13.71) . we will use the following 

Lemma 3.1 z(x,t) — > fn locally uniformly in Q as t — > +00. 

Indeed, from Lemma I3TTI since K > |m|, it easily follows that w(x,t) — > m locally 
uniformly in Q as t — > +00, and we get (13.71) with C = m. 

Proof of Lemma 13.11 We split it into several parts. 

1. Let iefi and r > such that B(x, r) C f2. From Proposition 13.11 it follows that 
there exists a sequence t n — > +00 as n —>■ +00 and a function z such that 

z(x,t n ) — > z{x) uniformly on B(x,r). 
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We define 

z n (x, t) := z(x, t + t n ) for all (x, t) G Q, x (— t n , +00) 

and notice that z n is a viscosity solution of (13.10!) on Q x (— t n , +00). From the uniform 
boundedness of z, we derive the one of z n and the half-relaxed limits method implies 
that 

z(x,t) — limsup z n (y,s) 

n — » + oo 

is a viscosity subsolution of the generalized Dirichlet problem 

Cz(x,t) <0 in Q x (— 00, +00) , > 

min{£z(x, t), z(x, t) + K} <0 on d£l x (— 00, +00) ^ ' ' 

2. We claim that max z(x, t) = m. Indeed, on one hand, for all n e N, there exists 
x n G f2 such that 

z(x n , t + t n ) = max z(x, t + t n ) — max ,2 n (x, t) . 

From the compactness of Q, we have, up to subsequence, x n — > x for some igSI. It 
follows that 

m = lim m{t + t n ) = lim [max z(x, t + t n )] 

= lim z n (x n ,t) 

n— >+ix) 

< 5(5, t) < max t). 

On the other hand, for all a; G fi, we have 

z(x,t) = limsup z n (y, s) < limsup m(s + t n ) = m 

(y,s)-*(x,t) (s,n)~*(t,+oo) 

and the claim is proved. 

3. Let to > be fixed; as a consequence of step 2, there exists a point xq G VL such 
that 

max z(x, t ) = z(xo,t ) = m. (3.15) 
xen 

To end this proof, it is enough to prove that (for some t > 0) xo lies inside Q. 

Indeed, if Xq G Q, then the Parabolic Strong Maximum Principle and (13.151) imply 

z = m in x (— 00, t ] (3.16) 
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Therefore, by taking (I3.16P into account, for any x G B(x,r), it would follow that 
m — z(x, 0) = limsup z n (y,s) 



TL — * + OC 



limsup z(y,s + t n ) 

(y,s)-*(x,0) 



n — > + oo 



lim sup [z(y, s + t n ) - z(y, t n )] + lim sup z(y, t Tl 

(y,a)-+(x,0) 2/->a= 



Tl — > -\-OC 



71 — >- + 00 



< lim sup 

(j/,a)-»(x,0) 



l^'-)llo!S(M-)l'l f 



+ z(aj) = z(x) 



with the last inequality following from Proposition !3.H -(ii). Finally, since z is inde- 
pendent of the sequence (t n ) n , we have obtained 

z(x, t) — > m uniformly on B(x, r) as t — ► +00, 

for any x £ Q and r > such that B(x,r) C fi, thus proving the statement of 
Lemma 13.11 

4. Now, we are going to prove that, for some to > 0, z(x, to) has a maximum point 
xo in Q. For that purpose, let us assume by contradiction that 

xq E d£l and z < fn in f2 x (0, to]. (3-17) 

We pick some < S < to and argue in the subset VI s := {x E £1 : < d(x) < 5} by 
introducing the function x defined by 

X(x, t) :=m + k(e~ d{x) - 1) - k(t - t ) for all (x, t) et¥ x [t - 5, t ] 

where 5 and k are to be chosen later in such a way to obtain 

z < X in ^Fx [t -5, to]. (3.18) 

(i) We start by proving that x satisfies 

Cx(x, t) > in tt 5 x (t - 5, <„]. 

Computing we have 

£ x ( x ,t) = -fcd(x) + fce" d(:c) [d(a;)Arf(a;) - - ■ Dd(x)]. 

By using (I3.13p . we deduce that 

C x (x, t) = ke~ d(x) [-d{x) e d{x) + d(x)(Ad(x) - 1) + m/(m - 1) + o*(l)] . 
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By the regularity of dfl, Ad is bounded in fl s and since d{x) < 5, we can choose 5 > 
small enough such that 

-d{x) e d(x) + d(x)(Ad(x) - 1) + m/(m - 1) + o s {\) > , 

for all x such that d(x) < 5. Hence we obtain 

Cx > on tt 5 x (t - 5, t ] . (3.19) 

(ii) For all (x, t) € <9f2 x (to — 5, to], we have x(^, t) =m — k(t — t ) > m and knowing 
that i < m in x 1, we easily conclude that 

z<x on dQ x (t - 5,t ]. (3.20) 

(hi) We set ^ = {i 6 fl : ct(x) = 5}. For any (x, t) G Tj x (t — 5, to], we have 

X(x, t)=m + k{e~ s - 1) - fc(t - t ) > m + k(e' 5 - 1). 

Moreover, since z(x, i) < m in Q x (0, to], we have 5(x, t) < m on x [to — 5, t ]. We 
use the upper semi-continuity of z to define 

r)(8) := min (m — 5) 

r a x[i -<5,t ] 

and we find that z(x,t) < m — on Ts x [t — 5, to]. By choosing > such that 
k(e~ s - 1) > -77(6), that is 

0<fc< (3.21) 

one gets 

5< X onV s x(t -S,to]. (3.22) 

(iv) For all x G we have 

X(x, t - 5) = m + A;(e" d{a;) - 1) + Atf > m + £;(e~ 5 - 1 + 5) > m 
hence we deduce 

z < x on fF x {t - 5}. (3.23) 

5. Now we use (jSJjj) in the interval (t - S,t ), together with (15351) . fl3T20|) . (ET22]) 
and ( 13.231) , so that applying the comparison principle on fl s x [to — 5, t ] we finally 
conclude that ( 13. 18ft holds. 
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Since z(x ,t ) = m = x( x o,to), h follows from (13.181) that z — x achieves its 
global maximum on VI s x [t — 5, t ] at (xo, t ). By using the definition of the viscosity 
subsolution z of (13.141) . it necessarily follows that 

mm{£x(x Q ,t ),m + K} < 0. (3.24) 

But Cx(xo, t ) = > and m + K > 0, hence we reach a contradiction with (I3.24p 
and we conclude that (13.171) does not hold. Hence, we have that (I3.15P holds for some 
t > and x G Q, and by Step 3 we deduce that (13.161) holds and the end of the 
proof follows. 

□ 



4 The Non- Convergence Cases: c = or c > and 

1< m < 3/2 

The main result of this section is 



Theorem 4.1 Let f2 be star-shaped with respect to a point xq G M. N . Assumqj that 
there exist constants 5 > and fiQ > such that we have, for any r G [1 — /io, 1)-' 

f(x)+c + lf{X0 + r{ *- X ^- nx)l <-5 on H. (4.1) 
(1-r) 

If c = and 1 < m < 2, or if c > and 1 < m < | ; i/ien 

u(x,t) + ct — > — oo as t —> +oo locally uniformly in Q. (4-2) 

More precisely, there exist a solution ip of U.8\) and continuous functions r(t), H(t) 
such that 

2 777, 

w(x,t) + ct < r(t)^iy?(xo + r(t)(x-x )) - (4.3) 
where r(t) | 1 and i?(t) ~~ * +°° as t — > +oo mf/i t/je following rate 



(a) Ifc> t/jen 



H(t) = 0{r^) ifl<m< 3 



2 ' 



H{t) = O (log t) i/m 



3 
2 ' 



(4.4) 



2 this assumption is for instance verified when / 6 W /1,oc (f2) and 

/(z)+c+||I}/|| 00 diam(ft) < -<5 
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(b) If c = then 



H(t) = 0{f- m ) ifl<m<2, 
H(t) = 0(\ogt) ifm = 2. 



(4.5) 



Let us recall that the ergodic constant c depends itself on /, so that the as- 
sumptions made in the above theorem are not obviously checked. However, such 
assumptions can actually happen to be true. To fix the ideas, consider the following 
Example which provides a specific case where condition (14.11) holds. 

Remark 4.1 Assumption ( f^.ip is always verified when f is a constant. 

We first consider the case that f = and denote by Cq the corresponding ergodic 
constant. One can easily deduce that c < 0. This is a consequence of the charac- 
terization \2. 1 6}) in Proposition \2.2\ indeed the constant functions are subsolutions 
of Ii2.15\) with A = while Proposition \2.2\ yields that the infimum in 112.16]) is not 
attained. Consider now the case that f — f is a constant, possibly different from 
zero. Of course the corresponding ergodic constant now is c = c — fo- For every 
value of f , (Jjjty is clearly verified with any 5 < |c |. 

Next, it is not difficult to construct some function f which is not constant and 
such that Iji4-1\ ) is verified. In particular, a small perturbation (in Lipschitz norm) of 
the constant fo still verifies fl^.lp ; indeed it is easy to check that the ergodic constant 
c depends continuously on f (with respect to perturbations in the sup-norm) as a 
consequence of formula K2.16) . 

Proof of Theorem 14.11 Without loss of generality, assume that xo = 0. Moreover 
observe that we can always replace 5 with a smaller value in (14. ip without loss of 
generality. 

2 rn 

In order to prove that (14. 3p holds, it is sufficient to prove that rlt)™^ p(r(t)x) — 
H(t) is a supersolution of E(Q, f + c, g + ct,uo), since then the estimate ( 14.31) will 
follow by an application of the comparison result. Let then p be a solution of (11.81) . 
we define the function <p r on Q x [0, +oo) as follows: 

2 rn 

ip r (x,t) = r(t) m - 1 ip(r(t)x) 

where r(t) will be chosen below in a way that < r(t) < 1 and r(t) j 1 as t — > +oo. 
In particular, we fix 

r(0) = l-/i 

so that r(t) G [1 — /i , 1) and (14.11) may be applied. Moreover, by the uniqueness 
result for the ergodic problem, we can represent ip as ip = tpo + L with p Q being the 
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unique solution of (jl.8p such that min^<£> = and L a constant to be fixed later. In 
particular, note that ip > L since (po > on Q. 

1. Our first step is in Q. It is easy to see that 

-A<p r + \D<p r \ m = {-Aif + \D<p\ m ) 
hence, with the notation A(z) = zt — Az + |-D,z| m , we get 

2 — m <p(r(t)x) 



A(<p r ) = r{t) — (f{r{t)x) + c) + r(t)^r'(t) 



m — 1 r{t) 



+ D(p(r(t)x) • x 1 . 



(4.6) 

Using the result (12.41) . we know that = o ^g^yj as x — > <9f2, where r and z/ 

are tangential and normal vectors. Then, since f2 is star-shaped with respect to 
and since Dtp ■ v > 0, we deduce the existence of o > such that < a implies 
Dip{y) • y > 0. Hence, since is smooth inside f2, we have 

Dip(y) ■ y > X{d(y)>*} > -C 
for some constant C a > 0. We deduce that 



Dcp(r(t)x) ■ x > 

In particular we have, for L large enough, 
2 — in (p(r(t)x) 



r(t) ' 



m 



+ D<p(r(t)x) ■ x > 



1 



2 — m 



L-CA > 



1 r(t) ' r(t) \m — 1 

Going back to ( 14.61) . we have, using that r(t) < 1, 

A(<p r ) > r(t)^(f(r(t)x)+c) 

= r(t)^(f(r(t)x) - f(x)) + (r(t)^r - l)(/(x) + c) + /(z) + c 



> 



r(t)) l/(^(^) /(x))| + ^ _ + c) + + c 



> -(l-r(t)^r) 



(1 -»•(*)) 
|/(r(f)s) -/(*))! 



(l-r(t)) 

Now we use (14.11) and we finally obtain 

A(<p r )>f(x) + c + 5(l-r(t) 



f{x) + c 



+ f(x) + c. 



(4.7) 
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Setting 

H(t) = (l-r(s)^) ds, 

it follows that 

A(<p r (x, t)-5 H(t)) > f(x) + c = A(u(x, t) + ct). 

2. Now, we consider the boundary dfl and we turn to the choice of r(t), which 
depends on the different values of m and c. 

(i) When c > and 1 < m < |. First observe that from (12. ip . when x € dQ, we 
have that 

2 — m 

<p(r(t)x) > K*(l -r(t))~^ + L (4.8) 
for some positive constant K*. Then we choose r such that 

m— 1 

(..- 2-m \ 2-m 
C + ^ + fft.- J for a.l t >0 

where A > will be determined later on. Note that r(t) is increasing, r(0) = 
1 — /i , and r(t) "f 1 as £ — > +oo. On the boundary, we get 



2 — m 



¥>rOM) > (l-/u )~ ( (c + A)t + L + K*n m ~ 1 
hence, up to choosing A sufficiently large (only depending on c and /io), we have 

2- / 2 — m \ 

yvfot) > 2c* + (1 -//<,) ^ Ll^o" 1 "' • ( 4 - 9 ) 



Next, since 



H(t) = J (l- r(s)^) ds<Kt 



for some if depending on c, if*, /io, A, and since the boundary datum g is 
bounded on dQ, we have from (14. 9p 

(p r (x,t) > g(x) + ct + 6H(t) 

up to taking in (14. ip some 5 eventually smaller and choosing L large enough. 
We conclude that 

(p r (x, t) — H(t) > u(x, t) + ct for all x G <9f2, t > , 

where = 5 H(t). 
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(ii) When c = and 1 < m < 2, we set 

l-r(t) = (/io^+t)- (m - 1} 

As before, r(t) is increasing, r(0) = 1 — ^0 and r(t) j 1. Using (14.81) we have, 
for every x G <9Q, 

^ r (a;,t) >r(0)%%<p{r(t)x) > (l-// )^ (V^ ^ + i) 2 ~ m + L 

while 

= f f 1 - r(s)^) ds<K f\l- r(s)) ds = ^ + tf~ m . 

Jo ^ ' Jo 2-m 

Therefore, choosing 5 eventually smaller in (14.11) . and then choosing L suffi- 
ciently large, we obtain again 

<p r (x, t) > g(x) + 5H(t) for all x G dVt, t > , (4.10) 

which means that <p r (x, t) > u(x, t) + H(t) on the lateral boundary with H(t) = 
5H(t). 

(iii) When c = and m = 2. Here recall that we have 

<p(r(t)x) > - log(l - r(t)) -K* + L 
for some positive constant K*. We set now 

l-r(t) = — L-j, 

hence 

¥V(ar, t) > log(t + /i : ) -K* + L. 

Since 

F(t) = J (l- r(s)^t) ds<K (1 - r(s)) ds = K \og{t + ^ x ) , 
we get again (14.101) by choosing 5 sufficiently small and L large enough. 

3. Finally, at £ = 0, we have 

2 — m 2 — m 2 — m 

<f r (x, 0) = r(0)^^(r(0)a;) = (1 - /i )^<^((l - /U )x) > L(l - // )^ 
hence for L large we obtain <p r (x, 0) > Mo(^) for all x G O. 

4. We conclude from the above comparison that <p r (x, t) > u(x, t) + ct + Now, 
according to the value of c and m, the choice of r(t) in (i)-(iii) gives the rates (14. 4p 
or (14. 5 p claimed for H. □ 
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5 Appendix 



In this section, we will give the proofs of Theorem 13.31 and Proposition 13. H - (i) . 
Proof of Theorem 13.31 

1. Let d(x) be the signed distance function, which is negative when x (jL Q. Let us 
fix 5 > such that d(x) is C 2 in {x G M. N : \d(x)\ < 25 } and set 



So 
2 



d(x) = x{d{x)) , 

where x{ s ) is a smooth, nondecreasing function such that x( s ) = s f° r < |s| < 
and x( s ) 1S constant for \s\ > 5q. Without loss of generality, we may have x'i s ) < 1 
for every s. Consider now the vector field 

n k (x) = - Dd(y) p k (x - y) dy 
Jr n 

where p k is a standard mollifying kernel (supported in the ball Bi(0)). Recalling the 

definition of d, the field n^(x) is supported in a neighborhood of dQ of radius 5q + 4, 
and we have G C°° for k large. Moreover, using the properties of d(x) and in 
particular that d G C 2 , we have 

KI<1, \Dn k \ < H^dlloo, \D 2 n k \ < . (5.1) 

Clearly is an approximation of the normal vector u(x) = —Dd(x), in particular 

\\D 2 d\\ 

\n k (x) + .Drf(x) I < " , "°° . (5.2) 

/c 

Then we consider the function 

v(x, t) = j(t)ip(x — n(t) n k (x)) where <f(x) = ifo(x) — L. 

Here (p is the unique solution of (11.81) such that miny9 = 0, and L is an additive 
constant to be chosen, whereas j(t) and pit) are positive functions, with values in 
(0, 1), that will be fixed later in a way that ^(t) j 1 and p(t) | as t — ■> +oo. Observe 
that since 

d(x — /i(t) nfc(x)) = <i(x) — pDd(x) ■ n k (x) + 0(p 2 ), 

by choosing k > 2||.D 2 <i|| 00 it follows that for any i6ll such that d(x) < y, we have, 
by using (15.21) and the definition of d, that: 

1 3 

d(x) + -p + 0(p 2 ) < d(x - p(t) n k {x)) < d(x) + ^ + 0(p 2 ). (5.3) 
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In the following, we fix k as above. Moreover, in order to have /j,(t) sufficiently small, 
it is enough to fix /x(0) small enough, since fi(t) is decreasing. To fix the ideas, we 
set (3 = yu(O), and we choose (3 small enough so that, thanks to (15.31) . we have the 
following: 

x — u(t) n k (x) ef2 Vi £ 0,t > 0, 

(5.4) 

d(x) + \/i(t) < d(x - /i(t) n k (x)) < d(x) + 2//(t) in tt s °/ 2 x (0, +oo). 

Note that this choice of (3 only depends on 5 and H-D^Hoo, in other words only on 
the domain Q, and, eventually, we are allowed to take a smaller value of (3 if needed 
later. 



2. Let us compute now the equation for v. Henceforth, we denote, for any function 
z, A(z) := Zt — Az + \Dz\ m , and we use the letter K to denote possibly different 
constants only depending on Q, m, f, c. 

Since Dv = "f(t)(I — }J,(t)Dn k (x))D<p(x— /i(t) n k (x)), using (15. ip and that 7, /i < 1, 
we get 

A(v) < --fAip + 7 m (1 + Kfi) m \D<p\ m 
+ 1 K^{\Dn k {x)\ \D 2 ip\ + \Dip\)+i(t)ip- ^{t)^ Dip ■ n k {x) 
where the argument of p is x — fj,(t)n k (x). Since p satisfies (11.81) . we deduce that 

A(v) < j(f(x - n(t)n k (x)) + c) + 7 ( 7 m - 1 (l + K/j) m - 1) \D<p\ m 
+ 7 Kfi (\Dn k (x)\ \D 2 p\ + \D<p\) + i(t)p - /i'(t) j(t)Dp ■ n k (x). ^ 

From the Lipschitz continuity of /, we have 

7 (/(x + fi(t)n k (x)) + c)< (f(x) + c) + (1 - 7 (t))||(/ + c)-\\ L «, {n) + fi(t)\\Df\\ L ^ {n) . 

Moreover, since p = <p — L, using (12.71) for p and that 2 — m < 1 we deduce that 
we have, in the whole range 1 < m < 2, 

p < K(l + \Dp\) - L VxGfi, 

for some constant K > 0. We also use (12.61) . which we can suppose to hold true in 
the support of n k (x) without loss of generality. Therefore, we obtain from (15.51) : 

A(v) < /(o;) + c+(l-7(t))||(/ + C )-|Uc« (n3 + / x(t)|| J D/||^ (n) 
+7 ( 7 m " 1 (l + Kfi) m -1 + K/i) \Dp\ m 
+j'(t)[K(l + \Dtp\) — L]+ 7 (f) K(\iS(t)\ + Kt))\Dp\ , 
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hence there exists a constant, still denoted by K, such that 

A(v) < /(x) + c+(l-7(t))||(/ + c)-|| £ co ( n)+M*)IP/IU-(n) 
+7 (7™" 1 -1 + Kfi) \D<p\ m + i{t) [K(l + \D<p\) - L] 

Here we take L > K and we choose j(t) such that 

7 (t) m - 1 = 1 - A/i(f) (5.6) 

for some A large enough. Without loss of generality, we can assume that fi(t) is small 
in a way that j(t) > (this amounts to ask (3 = /i(0) < v). Then we obtain 

A(v) < /( x ) +c +(l- 7 ( i ))||(/ + c )-|| ioo(n) + ^)||D/|| i oo (n) 

+ 7 /i (ff - \)\D<p\ m + i{t) K \Dtp\ + 7 (f) K{\pi\t)\ + //(f))|^| 

which yields, by applying Young's inequality, 

A(v) < /(x)+c+(l-7(*))||(/ + c)-|Ucc (n) +^)||I>/|Ucc ( n ) 

m 

+±-<yn(K-\)\D<p\ m + K>yfi ( ^ ) +K lf i ^ ! '' r| 



2 \ 7 M / V V 

hence, choosing A > K we get 

< /(x) + c +(l- 7 (t))||(/ + c)-|Uoc (n) +^11,0/11^(0) 

m m 

+A - 7 ,f2wy"- +/f7 ,f^V-. (5 . 7) 



IV- J \ v> 

3. Let us consider now the boundary. Since we have from (15.41) 

d(x — (j,(t)rik(x)) > T^(f) f° r an x ^ ^ ) 

using the asymptotic behavior (I2.1I) - (I2.2I) for (p we deduce that there exists a constant 
K* such that 

<p Q (x - n{t)n k {x)) < K*(fi(t)- a + 1) - L for all 
when a = ^7 > (i.e. m < 2), while when m = 2 we have 

m— 1 V /' 

(po(x — jj(t)rik(x)) < — log n(t) + If* — L for all x G <9f2 . 
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We are going to take later L sufficiently large. In particular, considering L > K* we 
deduce that for every x e dfl and t > 0: 



if 1< m < 2, v(x, t) < K*fx(t)~ a + i(0)(K* - L) 
if m = 2, v(ac, t) < - log /x(t) + 7(0) (iT - L) 



(5.f 



Moreover, for the initial condition we have 

v(x, 0) = 7(0)<^o(x - fi(0)n k {x)) - 7 (0)L . (5.9) 
Let us now distinguish the choice of fi(t) according to different situations: 
(i) If c > and 1 < m < 2, we set 

'ct + K*(3- aS 



K* 



where a = ^zj- The value of (3 = /i(0) has been already chosen, as explained 
before, in order that f!5.4p holds true and also 1 — \[3 > 0, this choice only 
depends on Q, f, m. Note that fi(t) is a decreasing function and satisfies 

K*^(t)- a = ct + K*p- a } 

hence ( 15.81) implies 

v(x,t) <ct + K*f3~ a + 7(0) (K* - L) for all x e dVt . 

Since 7(0) m_1 = 1 — A /i(0) = 1 — \(3 > 0, and since g is bounded in dfl, up to 
choosing L sufficiently large we will have 

v{x,t) < g(x) + ct for all x E dtt. (5.10) 

Similarly, from the boundedness of Uq on Q, up to choosing L large enough 
(again depending on /3), we obtain from (15. 9p 

v(x, 0) < Uo(x) for all x G Q. (5-H) 

Consider now the behavior of /J,(t), which implies that fi(t) — > 0, and fi'(t) = 
o(/i(t)) as t — » +00. By definition of 7 in (15.61) . we have that j(t) | 1, 1 — 7(t) = 
0(/x(t)) and 7 ; (i) = 0(|/x'(t)|) as t — > +cxd. In particular we have 

7 K^)"' %7 n^)""- A> (5 ' 12) 
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hence we deduce from (15. 7p 



A(v) < f(x) + c + Kfi(t) (1 + \\(f + c)-|| L <» ( n) + ||I>/IU~(n)) . (5.13) 
Therefore, if we set 

K = K(l + \\{f + c)-|U=o (n) + ||D/|| L <» (n) ) 

and we define 

w = u + ct + K / fj,(s)ds , 



we can conclude that > A(v). Moreover, from (15.101) and (15.111) we have 

that «<toon the parabolic boundary. By the standard comparison result, it 
follows that 

u(x, t) + ct> v(x, t) —K /J,(s)ds (5.14) 

Jo 

for all (x,t) G Q x [0, +oo). Since /i(t) = O as i — > +oo, the conclusion 

of (13 .4p follows according to the values of a, i.e. of m. Note in particular that 
/J,(t) G -^ 1 (0, oo) if and only if a < 1 which corresponds to m > |. 

If c > and m = 2, we set 

Mt) = /5e- ct , 

where /? = /x(0) is chosen as before. With this choice we have from (15.81) 

v(x,t) < ct -log/3 + 7(0)(iT - L) , 

and choosing L large enogh we deduce (15. lUj) . Of course (15.111) remains true 
as before. Finally, in this case we have jj!(t) = 0(fi(t)) as t — ► +oo, and then 
7'(t) = 0(/i(t)) as well; therefore (15.121) still holds true. Thus we obtain again 
(15.131) and we conclude as before the inequality (I5.14p . Being // integrable in 
(0, +oo), this implies (13.41) for the case m = 2. 

If c = and 1 < m < 2, we set as 

where A will be fixed later. The initial condition /i(0) = /3 is fixed as before. 
We have 

,/(*) = - A^ff - < A/x(*) 
^ v y /3~( 1+a ) + At ^ w 
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hence, using also the definition of 7, we obtain (15.121) with a constant depending 
on A. We deduce then from (15.71) (recall that here c = 0) 

A(v) < f(x)+Kti(t)+K At i(t), (5.15) 

where 

K = K (1 + ||(/)-|Uoo (n) + Ip/Hoo) 

is a constant only depending on Q, m, f. In particular we obtain that the 
function ^ 

w = u + (K + K A ) / /i(s)ds (5.16) 
Jo 

satisfies A{w) > A(v). Now we choose A = , in order to have 

K*fi(t)~ a = K* (p- {1+a) + At)^ 

= K *A^— [ fi(s) ds + K*/3~ a 
1 + ot Jo 

= K I fx(s)ds + K*P~ a 
Jo 

< w(x,t) -u(x,t) + K*(3~ a . 

Therefore from (15. 8p we deduce that on the boundary 

v(x, t) < w(x, t) - g(x) + K*/3~ a + 7(0)(iT - L) for all x G dVl. 

Choosing L large enough we conclude that v(x,t) < w(x,t) for every x E dQ 
and if: > 0; as before, we also have (15. lip , hence by comparison we conclude that 

u{x,t) > v{x,t) - (K + Ka) [ fi{s)ds 

Jo 



Since J Q * /i(s) ds — O f^ 1 " 4 "™) as t — > +00, we obtain (13.51) by definition of 
If c = and m = 2 we set 

1 



a. 



At + p- 1 
Since we have 
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we obtain (15.151) as before, hence A(w) > A(v) where w is defined in ( 15.161) . 
Since t 

— \ogfi = A / fi(s) ds — log (3 
Jo 

we choose A = K and by using (15.81) we obtain on the lateral boundary: 

v(x,t) < K [ fi(s) ds- log + 7(0) (iT - L) 
Jo 

< w(x, t)-g{x)- log (3 + j{0){K* - L) . 

We conclude that d<woii the boundary up to choosing L large enough and 
then, by comparison, we deduce that v(x,t) < w(x,t) in Q x (0,T), which 
implies ( 13.51) for the case m = 2. 

□ 

Now, we turn to the proof of the interior Holder estimates on v with respect to the 
a;-variable uniformly in t 6 (0, +oo) which is based on an idea introduced by Ishii and 
Lions [12]. This idea has been already used for instance in Barles pQ and Barles and 
Souganidis [7] to show gradient estimates of viscosity solutions to quasilinear elliptic 
ans parabolic PDE with Lipschitz initial conditions, by Barles and Da Lio [5] to prove 
local Holder estimates up to the boundary of bounded solutions to fully non linear 
elliptic PDE with Neumann boundary conditions and by Da Lio [10] to obtain C 0,u - 
estimates for viscosity solutions of parabolic equations with nonlinear Neumann-type 
boundary conditions. 

Proof of Proposition !3.H -(i). We fix v 6 (0, 1), 5 > 0, rj > and x G such that 
d(xo) > 5. We are going to show that there exists a suitable constant C depending 
on z/, 5, i], M := sup KsX ^/ 2 ,+oo) M an d the data of the problem such that, for all 

y G B(x , 5/2) and t > 77, we have 

v(x , t ) - v(y, t ) < C\x - y\ v . (5.17) 

The property (15.171) clearly implies the C^-estimates of v(-,t ) in K$. Indeed, for all 
x, y G Ks, if \x — y\ > 5/2 then 

2 u+l M 

\v(x,t ) -v(y,t )\ < — \x-y\ v . 

But when \x — y\ < 5/2, we can therefore apply (15.1 Tf) with xo = x and y G B(x,5/2) 
or with xo = y and x G B(y, 5/2) and finally obtain the desired estimates. 
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To prove (15.171) . we consider the function (x,y,t) i— »■ $(x, y, t) defined on Q x Q x 
(0, +00) as follows 

®(x,y,t) = v(x,t) - v(y,t) - <p(\x - y\) - L(\x - x \ 2 + \t - t \ 2 ) 

where ip(t) = Ct u . The constants L > and C > will be chosen in such a way that 
$ is a non-positive function. 

We first choose L and C in order to have 

$<0 on d[B(x , 8/2) x B(x , 5/2) x (77/2, +00 



This leads to the constraints 



L (Q >2M, L(|) 2 >2Mandc(£) > 2M . 

With these choices of C and L, it is easily checked that §(x,y,t) < for any x,y,t 
such that \x — y\ > 5/4 or \x — x \ > 5/4 or \t — t Q \ > r]/2, and putting together these 
properties, we clearly have the desired property. 

Next, L being fixed as above, we argue by contradiction, assuming that, for any 
constant C (satisfying the above constraints), we have 

M c ,l ■= _ _max $(x,y,t)>0. (5.18) 

B(i ,ii/2)xB(io,«/2)x[>)/2,+oo) 

Let (x, y,i) e B(x , 5/2) x B{xq, 5/2) x [r//2, +00) be a maximum point of we have 
dropped the dependence of x,y and t in C for sake of simplicity. By using (15 . 1 8j) . it 
is clear that x 7^ y, otherwise, we would have M.c,l < 0. With the choices of C and 
L we made above, it is obvious that (x,y,t) e B(x ,5/2) x B(x ,5/2) x (r//2,+oo). 
From (I5.18p . we get 

C\x - y\ v + L(\x - x | 2 + |t-to| 2 ) < v(x,t) -v{y,i) (5.19) 

which yields 

C\x-y\ u <2M and L(|x - x \ 2 + |f - t \ 2 ) < 2M (5.20) 

and it follows, in particular, that \x — y\ — * as C — > +00 and we recall that we may 
assume without loss of generality that \x — y\ > for C large enough. 
We define by 

<j>(x-y) = ip(\x-y\) (5.21) 
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and use the arguments of [121 Proposition IV. 1] to prove the existence of two N x N 
symmetric matrices Bi and B 2 , and a, b G M such that 

B 1 -2LI \ / D 2 <\>{x-y) -D 2 <P(x-y)\ 

-B 2 J ~ \ -D 2 (f)(x - y) D 2 cf){x - y) ) ' l ° } 

a-b-2L(t-t ) > and 

a - Tt(B x ) + \p\ m - f(x) - c < 

6-Tr(5 2 ) + |g| m -/(y)-c>0 1 J 

where 

p = D<f)(x — y) + 2L(x — x ) and g = D<fi(x — y). 
For all £, C e M^, we rewrite (I5"^2D as 

(£x£, - (S2C, < (^ 2 0(^ - - 0, £ - + 2L|£| 2 . (5.24) 

Let (ej)i<j<7v_i be a familly of (iV — 1) vectors in ~R N such that (ei, e 2 , ejv-i, 1%) is 
an orthonormal basis of M. N . Plugging successively £ = C = e « f° r all i = 1, iV — 1 
and £ = — £ = A in (I5.24p and by adding all the inequalities obtained, we obtain: 

2V-1 

Tr^-Bi) = y2((Bi- B 2 )e t , ei ) + {{B t - B 2 )^,^} 

< 2iVL + 4 (D2 ^7/ )g ' g) . (5.25) 

Going back to the form of 0, we set x( z ) = \ z \ an d obtain 

(D 2 (f)(x - y)q, q) = (f'{D 2 x(x - y)q, q) + (p"((Dx(x - y) ® - y))q, q). (5.26) 

But, knowing that D\{z) = A, it follows that |.Dx(z)| 2 = 1 which yields by differ- 
entiation 2D 2 x{z)Dx{z) = and we obtain D 2 x{x — y)q = by taking z = x — y. 
Therefore, we use (15. 26ft and find that (15. 25ft becomes 

Tr(B 1 - B 2 ) < 2NL + ^"iii^Ml = 2 NL + Ay" 

\q\ 

Finally, we obtain 

Tr(B l - B 2 ) < 2LN + ACv{v -l)\x- y\ u ~ 2 . (5.27) 
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By subtracting the two inequalities in (I5.23P and using (15.271) . we have 

2L{i - 1 ) + 4CV(1 -v)\x- yr 2 < \q\ m - \p\ m + f{x) - f{y) + 2LN . 
Then the convexity of p i— > \p\ m yields \p\ m > \q\ m + m\q\ m ~ 2 q ■ (p — q) and therefore 
\q\ m ~ \p\ m < ml^r^b - q\ < mLSlql™- 1 , 

since \x — Xq\ < 8/2. 

Plugging this estimate in the above inequality gives 

2L(t-t ) + 4CV(1 -u)\x- y\ v ~ 2 < mU\q\ m ~ x + + 2LN , 

and therefore 

2L{t - t ) + 4CV(1 - - y\ u ~ 2 < mU[vC\x - y^)^ 1 + 211/11^ + 2LN , 

Hence, using that 2L\t-t \ < 2L 1 / 2 (2M) 1 / 2 , if K := 2L 1 / 2 (2M) 1 / 2 + 211/H^ + 2LN , 
we finally have 

4z/(l - u) < K ^~f} 2 - + mlfo™- 1 ^ 2 !! - y|("-i)("»-2)+i . 

But m - 2 < and {y - l)(m - 2) + 1 > 1 since v - 1 < as well; from (J520D, 
we know that |x — y\ — > when C tends to infinity. Therefore it is clear that this 
inequality cannot hold for C large enough, thus contradicting (I5.18p . Then, if C is 
large enough, the estimate (I5.17P holds true. Examining this proof, we find that C 
depends on u, 5, r], M := sup KsX ^/ 2 ,+oo) M> ll/lloc an d the constants m, N. 
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